Partial order from disorder in a classical pyrochlore antiferromagnet 
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We investigate theoretically the phase diagram of a classical Heisenberg antiferromagnet on the 
pyrochlore lattice perturbed by a weak second-neighbor interaction J2. The huge ground state 
degeneracy of the nearest-neighbor Heisenberg spins is lifted by J2 and a magnetically ordered 
ground state sets in upon approaching zero temperature. We have found a new, partially ordered 
phase with collinear spins at finite temperatures for a ferromagnetic J2. In addition to a large 
nematic order parameter, this intermediate phase also exhibits a layered structure and a bond order 
that breaks the sublattice symmetry. Thermodynamic phase boundaries separating it from the fully 
disordered and magnetically ordered states scale as 1.87 J2S 2 and 0.26J25' 2 in the limit of small J2. 
The phase transitions are discontinuous. We analytically examine the local stability of the collinear 
state and obtain a boundary T ~ J%j J\ in agreement with Monte Carlo simulations. 



I. INTRODUCTION 

Magnets with geometrical frustration^, have received 
much attention as models of strongly interacting elec- 
tronic systems with unusual ground states, thermody- 
namic phases, and excitations. The hallmark of strong 
frustration is a conspicuously large degeneracy of the 
classical ground state: essentially, a finite fraction of the 
degrees of freedom remains unconstrained to the lowest 
temperatures. For discrete spins, this manifests itself in 
the number of ground states scaling exponentially with 
the system volume and thus giving rise to a nonzero en- 
tropy density at absolute zero temperature. Well-known 
examples of that are the Ising antiferromagnet on the 
triangular lattice^ and spin ice£ For continuous spins 
— most saliently for the Heisenberg antiferromagnet on 
the pyrochlore lattice — the classical ground states form 
a manifold whose dimension is proportional to the sys- 
tem volume^ In that particular case, the classical model 
exhibits strong short-range spin correlations but fails to 
exhibit any form of conventional magnetic order down to 
the lowest temperatures accessible in Monte Carlo sim- 
ulations. The strong correlation between the local mo- 
tions of spins in this liquid-like phase manifests itself as 
an emergent gauge structure in the low-temperature limit 
and results in a dipolar form of the asymptotic spin cor- 
relations at large separations^ 

At the same time, the large degeneracy of the ground 
state makes this system susceptible to all kinds of pertur- 
bations, which certainly exist in real compounds. For in- 
stance, the spin-lattice coupling, arising from the depen- 
dence of exchange strength on the atomic displacements^ 
lifts the degeneracy through a spin analog of the 
Jahn- Teller effect^ observed in spinels ZnC^O/t 1 — and 
CdCr 2 4 .-^ 

This naturally leads one to ponder the following ques- 
tions. Can the interplay of a weak perturbation with 
strong frustration lead to interesting ordered phases? 
Are there any (intermediate) partially ordered phases? 
What is the nature of the phase transitions between such 




FIG. 1: Second and third-neighbor pairs on the pyrochlore 
lattice. Since exchange paths giving rise to J3 and J3 are 
inequivalent, the two couplings may be different. Numbers 
from to 3 label the four fee sublattices. 



phases? In this paper we discuss these questions in the 
context of a classical Heisenberg antiferromagnet on the 
pyrochlore lattice with interactions going beyond near- 
est neighbors. Following previous work by Reimers et 
alH and by Tsuneishi et ai.p2 we consider the classi- 
cal Heisenberg antiferromagnet on the pyrochlore lattice 
with the Hamiltonian 
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where (ij) and ((ij)) indicate pairs of first and second 
neighbors, respectively. Given the short-range nature of 
exchange forces, we work in the limit J% -C J\. It is rea- 
sonable to expect that the influence of J2 becomes notice- 
able only at low temperatures of order J2S 2 , when the 
system is already in the strongly correlated paramagnetic 
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FIG. 2: Phase diagram of the model with antiferromagnetic 
first and weak second-neighbor exchange of either sign on the 
pyrochlore lattice. Open circles are numerically determined 
locations of thermodynamic phase transitions (all first order) ; 
filled circles denote the stability boundary of the collinear 
phase. Solid lines are interpolated phase boundaries; the 
dashed line is a boundary of local stability of the collinear 
phase. The wavenumber of the incommensurate magnetic 
phase is h ~ 3/4. 



state, in which it is constrained to fluctuate around the 
ground states of the nearest neighbor exchange. Using a 
combination of Monte Carlo simulations and analytical 
arguments, we have mapped out the phase diagram in 
the J2-T plane shown in Fig. [2] 

Antiferromagnetic second-neighbor exchange, J2 > 0, 
significantly reduces the frustration by selecting states 
in which spins within any of the four fee sublattices, 
comprising the pyrochlore lattice, are parallel to one an- 
other. We find a collinear ly ordered phase of the type 
(So) = —(Si) = — (S2) = (S3), where the subscripts enu- 
merate the fee sublattices [Fig. Q] and Fig. (3Ja)] . The 
transition between the paramagnetic and antiferromag- 
netic phases is discontinuous. 

Ferromagnetic second-neighbor exchange, J2 < 0, 
leaves the system strongly frustrated. A mean-field cal- 
culation by Reimers et al^ predicted a ground state 
with incommensurate magnetic order. While Tsuneishi 
et al^ indeed observed Bragg peaks in the spin struc- 
ture factor obtained through a Monte Carlo simulation 
for J2 = — O.lJi, they also noted that the spins remained 
dynamic, failing to freeze. We show that the observed 
locations of the Bragg peaks are compatible with the re- 
sults of Reimers et ai, so that the low-temperature phase 
is most likely magnetically ordered. 

The main focus of our paper is a peculiar partially or- 
dered phase sandwiched between the paramagnet and the 



magnetically ordered state for weak enough ferromag- 
netic J2, namely — 0.09Ji < J2 < 0. In the intermedi- 
ate phase, the spins display collinear order; furthermore, 
they exhibit magnetic order within a thin {100} layer but 
no order across different layers. The partial order can be 
characterized by a combination of a director n specify- 
ing a global spin axis, a Potts (Z3) variable q = (100), 
(010), or (001) specifying the direction of the layers, and 
an Ising {Z%) variable o~ n for each layer identifying one 
of the two possible spin orientations within a layer. The 
order is partial in the sense that the Ising variables {<r n } 
randomly pick values of +1 and —1 with no discernible 
correlations between adjacent layers. The partially or- 
dered state is bounded by first-order transitions on both 
the high and low-temperature sides. 

Similar partial order has been previously found in 
& 1/S treatment of the Heisenberg antiferromagnet on 
the checkerboard lattice, also known as the square lat- 
tice with crossings, a two-dimensional analog of the 
pyrochlore.— In both systems, the distinct layered states 
are not related to one another by a symmetry of the 
Hamiltonian and simply arise as different local minima 
of the free energy. Free-energy barriers separating them 
may be large enough in practice for the system not to 
be ergodic and instead to remain in one of these minima 
forever. 

Since the energy of the partially ordered collinear state 
is greater than that of the low-temperature multiple-q 
magnetic order, entropic selection plays a crucial role 
in the stabilization of the intermediate phase. This 
is consistent with the general observation that states 
with collinear spins tend to have softer thermal fluctu- 
ations and therefore have a lower free energy at finite 
temperatures^^ A similar collinear phase has been re- 
ported in the Monte Carlo study of a J- J' model which 
interpolates between the pyrochlore and the fee lattices^ 

While we have focused on the role of second-neighbor 
exchange J2 in the formation of magnetic order on the 
pyrochlore lattice, our results also shed light on the role 
of third-neighbor interactions J3 (see Fig. Q}. In view 
of strong correlations between nearest-neighbor spins de- 
veloping at temperatures well below J\S 2 , the properties 
of the system depend not on J2 and J3 separately but 
on their linear combination J 2 — J-$- Indeed, the rela- 
tive shift in energy for any pair of ground states of the 
nearest-neighbor exchange due to a small J3 is identical 
to the effect of a J2 of the same magnitude and opposite 
sign. Thus our findings should also be of relevance for the 
more general case of a pyrochlore antiferromagnet with 
small J2 and J3. 

The remainder of this paper is organized as follows. 
In Sec. |TT] we briefly discuss the nature of magnetically 
ordered phases at low temperatures for both signs of the 
second-neighbor coupling J 2 . Sec. [TTT] presents the main 
subject of this work, the partially ordered phase found at 
intermediate temperatures on the ferromagnetic side of 
J2 . Stability of the partially ordered state and its phase 
boundaries are examined in Sec. IIVI We conclude with a 
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FIG. 3: (a)A q = Neel order for model with an antifer- 
romagnetic J2 (ferromagnetic J3). The order parameter is 
one of the three staggered magnetization L3 = (So + Si — 
S2 — Ss)/4S.— (b) The phase transition between the param- 
agnetic and antiferromagnetic phases for J2 = O.OUi. The 
simulated system has 16 x 8 3 spins. The energy density 
e = (E — Eo)/6N s , where Eo = —N 3 Ji is the ground state 
energy of nearest-neighbor interactions. 

discussion of these results in Sec. [V] 



II. LOW-TEMPERATURE ORDERED PHASES 

Since the phase transitions shown in Fig. [2]are strongly 
discontinuous and occur at very low temperatures, the 
metastable states close to the coexisting region are rather 
long-lived. Conventional histogram methods with lo- 
cal Metropolis updates are ineffective in determining the 
critical points due to a large energy barrier separating the 
metastable state from the true ground state. Instead, 
we settled on using a method proposed by M. Creutz 
et al-r^ in which a mixed phase with the two coexist- 
ing states each occupying half the lattice is constructed 
first. By thermalizing the mixed phase at various tem- 
peratures, the critical point is determined when neither 
of the two states prevail the system during the relaxation 
process. Since the multiple-q magnetic order has an ex- 
tended unit cell with a period of about 4 cubic lattice 
constants, systems used in our mixed-phase simulations 
contain 8 3 cubic unit cells, with a total spin N — 16 x 8 3 . 



A. Antiferromagnetic J2: low frustration 

In the limit 3i <C J\ , magnetic ordering takes place at 
a temperature T c — 0(J2S 2 ). The nature of this ordering 
is best understood by appealing to the fact that a weak 
third-neighbor coupling J 3 <§; Ji (Fig. [J) selects among 
the nearest-neighbour ground states in the same way as 
a second-neighbor coupling J2 of the same strength and 
opposite sign, as explained in Appendix!^] (We here note 
in passing that, since the strength of coupling depends 
on the exchange paths and not the interatomic distance 
alone, sometimes J3 may be as big as Ji- For instance, ab 
initio calculations show that in CdC^CU J3 exceeds J2 



in magnitude. 18 ! 19 ) This insight is useful as the resulting 
ordered pattern can be understood in a more straightfor- 
ward way by analyzing the effect of J3. To see that, note 
that the pyrochlore lattice consists of four fee sublattices 
and that third neighbors on the pyrochlore lattice belong 
to the same fee sublattice (Fig.[T]). Thus a ferromagnetic 
exchange J3 < is not frustrated and will be absolutely 
minimized by a state where spins within the same fee 
sublattice are parallel to one another. 

A translationally invariant four-sublattice ground state 
was predicted for the pyrochlore antiferromagnet with a 
ferromagnetic J3 by Reimers et al*£ The same can be 
expected for an antiferromagnetic second-neighbor cou- 
pling J2 > 0. In both cases the energy of the further- 
neighbor exchange is minimized by a ferromagnetic order 
(Si) within the individual sublattices. Consequently any 
configuration satisfying 2j=o(^») = is a ground state 
at the mean-field level. Thermal fluctuations nonetheless 
favor those with collinear spins. 5 This is indeed what we 
obtained in the Monte Carlo simulations (Fig. [3]): a q = 
Neel state with an up-up-down-down spin configuration 
on every tetrahedron is found to be the ground state for 
an antiferromagnetic J2. This collinear magnetic state 
is separated by a discontinuous transition line from the 
high-temperature cooperative paramagnetic state. As 
shown in Fig. E^b), both the energy density e and the 
staggered magnetization L 3 = (So — Si — S 2 + S3) /4S 
show a clear jump at the transition temperature T c w 
3.2 J 2 S 2 . 



B. Ferromagnetic J2: high frustration 

The case of a ferromagnetic second-neighbor coupling, 
J2 < 0, is similar to that of J3 > 0. An antiferromagnetic 
coupling on an fee lattice is frustrated, so that this time 
one may expect a more complex magnetic order. Indeed, 
Reimers's mean-field calculation yields an incommensu- 
rate magnetic order with a wavevector q = 2w(h, h, 0) in 
the case of a ferromagnetic J2. 

We have performed Monte Carlo simulations on the py- 
rochlore lattice with periodic boundary conditions mea- 
suring 8 cubic unit cells in each direction. The simula- 
tions were done for J2 = —0.1 J\. They revealed a state 
with magnetic Bragg peaks at incommensurate lattice 
momenta near 27r{3/4, 3/4,0} and other equivalent po- 
sitions. Fig. [Ua) shows two inequivalent Bragg peaks, 
q « 2tt(3/4, 3/4,0) and -2tt(3/4, 3/4,0), the rest be- 
ing related to these two by a reciprocal lattice vector. 
Bragg peaks with comparable intensities are found at 
other wavevectors related to the above two by point- 
group symmetries. This multiple-q Neel order is consis- 
tent with the ground states of ([1]) in the spherical approx- 
imation, in which the local length constraints |Sj| = S 

are replaced by a global one Yli=i l^il 2 = NS 2 . Intro- 
ducing the Fourier transform Si = ^ q S m (q)e lq ri [the 
site index i = (m, r,), where m is the sublattice index], 
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FIG. 4: (a) The spin structure factor of the low-temperature 
ordered state at wavevectors q = 2ir(h,h,l). The state was 
obtained from a Monte-Carlo simulation with a system with 
16 x 8 3 spins for a ferromagnetic J2; the temperature was 
T — 0.2 \ J%\. (b) Minimum eigenvalue of the exchange matrix 
Jmn(q) at wavevectors q = 2ir(h,h,l) for J2 = — Ji/10. The 
satellite peaks at q ~ 2-7r(|, |, ±0.1) might be due to the 
finite-size effect for an incommensurate spin order. 

the exchange interaction (Q]) becomes 

AT 3 

U = ^T. E -Wq)S m (q)-S„(-q). (2) 

q m,n— 

The Fourier components S m (q) are subject only to a 
global constraint Yl m q |Sm(q)| 2 = $ 2 - The matrix 
J m „(q) is the Fourier transform of the exchange inter- 
action Jij — Jmn(^i — Tj). Its explicit form with interac- 
tions up to the fourth nearest neighbors can be found in 
Ref. El 

Expanding S m (q) = Ea^qm^q m terms of the 
eigenvectors U£ of the exchange matrix J m „(q) yields 
the energy as a function of the expansion coefficients <I? q : 

N 4 

q a— 1 

where A q is the corresponding eigenvalue of J mrl (q). 
With the normalization J2m=o l^q.ml 2 = 1' tne vectors 
3> q satisfy X^ q S a l* q l 2 = $ 2 - The ground state energy 
of © is thus E = NS 2 X min , where A min is the lowest 
eigenvalue A q . 




FIG. 5: The phase transition between the paramagnetic and 
antiferromagnetic phases for J2 = — O.lJi. The simulated 
system has a total of N = 16 x 8 3 spins. The normalized 
energy density e = (E — Eq)/6N, where Eo = —NJ1S 2 is the 
ground state energy of nearest-neighbor interactions. <j>M is 
the second moment of the magnetic order parameters. 

For the nearest-neighbor interaction only, the two low- 
est eigenvalues are q- independent, A q = A q = — Ji, re- 
flecting the degenerate nature of the magnetically or- 
dered ground state. This degeneracy is lifted by the intro- 
duction of J2 as discussed by Reimers et al^ 2 A contour 
plot of the lowest eigenvalue of the exchange matrix as 
a function of the wavevector q = 2n(h, h, I) for J 2 < 
is shown in Fig. 0Jb). ^ can be seen from Fig. [4] that 
the peaks of the spin structure factor appear at the same 
locations as the minima of exchange energy, namely at 12 
incommensurate wavevectors q* = 2ir{h*, h* , 0}, where 
h* » 3/4 depends weakly on the ratio J 2 j J\. For small 
JaM, h* =a + a 1 (J 2 /J 1 ) + O((J 2 /Ji) 2 ), where 

a = i arccos[(4\/3 - 9)/3] = 0.7427, 
44 

01 = = = 0.0336. (4) 

37T\/9654+ 5574V3 

The magnetic order is described by the order parame- 
ter composed of 12 vector amplitudes 3> q * A detailed 
characterization of this magnetic state is deferred to a 
future publication. Fig. [5] shows the temperature depen- 
dence of the energy density e and the magnitude of the 
order parameters <pM = l^q* I 2 - Both exhibit a clear 
jump at T c m 0.95 1 J 2 \S 2 , indicating a first-order transi- 
tion. This is also confirmed by a double-peak structure 
in the energy histogram at the transition temperature. 
Similar results were obtained for J 2 < —0.09 J\ where the 
magnetic phase is separated from the high-temperature 
spin liquid phase by a first-order phase transition as in- 
dicated in Fig. O 

III. PARTIALLY ORDERED PHASE 

As discussed in the Introduction, an intermediate 
phase with collinear spins exists at finite temperatures 
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FIG. 6: Transitions between (a) the paramagnetic and ne- 
matic phases, and (b) the nematic and Neel phases, for 
J-2 = — O.OUi. A parallel-tempering Monte Carlo method was 
employed to simulate a system with 16 x4 3 spins. The normal- 
ized energy density e = (E — Eq)/QN, where Eq = —NJiS 2 is 
the ground state energy of nearest-neighbor interactions. Q 
is the spin nematic order parameter. 



for a small ferromagnetic coupling J2 < 0. The ap- 
pearance of collinearity is not totally unexpected as it is 
well known that collinear states are in general favored by 
thermal fluctuations in magnets with frustrated exchange 
interactions.— The fact that the system remains frus- 
trated even in the presence of a ferromagnetic J2 makes 
the existence of the nematic phase possible. From an- 
other perspective, the classical nearest-neighbor Heisen- 
berg spins on the pyrochlore lattice evade the thermal 
selection only marginally^ The introduction of a ferro- 
magnetic Ji reduces the dimension of ground-state man- 
ifold, thus permitting thermal fluctuations to stabilize 
collinear states. 



A. Nematic order 

To demonstrate that spins indeed become collinear in 
the intermediate phase, we have obtained from Monte 
Carlo simulations the nematic order parameter Q de- 
fined as the largest eigenvalue of the traceless tensor 
Q\xv = (SfiS u /S 2 — Sfiu/S}^- where S 1 ^ represents Carte- 
sian components of a spin. It vanishes in a totally dis- 
ordered state and attains the maximal value of 2/3 for 
parallel spins. 

The thermodynamic behavior of the system with J2 = 




green 



FIG. 7: The three distinct collinear states of a tetrahedron. 
Frustrated bonds (with parallel spins) are shown as dashed 
lines. 



—0.01 J\ in the vicinity of the phase transitions is illus- 
trated in Fig. [6] The simulation was done on the py- 
rochlore lattice with periodic boundary conditions mea- 
suring 4 cubic unit cells in each direction, giving a total 
of N — 16 x 4 3 = 1024 spins. To improve the equilibra- 
tion process, we employed parallel temperin g 22 ' 23 with 
30 replicas. The energy density e and the nematic order 
parameter Q are shown as functions of temperature near 
T c i [paramagnet to partially ordered phase, Fig. [6] (a)] 
and T C 2 [partially ordered phase to antiferromagnet, Fig. 
[S] (b)]. The energy density shows a clear discontinuity at 
both transitions. Extrapolating the energy curve from 
the partially ordered phase to T = yields a density 
El = — 1/3 characteristic of a layered state to be dis- 
cussed below. Likewise, the order parameter Q extrapo- 
lates to the maximal attainable value of 2/3 characteristic 
of collinear spins. Below T C 2, the antiferromagnctic state 
seems to have a residual nematic order with Q ss 0.05, 
which may be intrinsic to the low-temperature ordered 
state, or a finite-size effect. 



B. Bond order 

Nematic order alone does not provide a full character- 
ization of this phase: four spins on a tetrahedron have 
three distinct collinear states not related to each other 
by a global rotation of the spins. They are labeled red, 
green, and blue in Fig. [Jj These states differ from one 
another by the location of frustrated bonds (ij) that in- 
volve parallel spins. Since the global direction of the 
spins is already captured by the nematic order parame- 
ter Q M „, further characterization can be made by using 
scalar quantities, such as bond variables /y = (Si • S^). 
At temperatures well below JiS* 2 only two (out of six) 
bond variables of a tetrahedron are independent: 9 



h 



h = 



foi + /23 + /02 + /l3 — 2/o3 — 2/i 



12 



foi + /23 — fo2 — f 



13 



(5) 



The vector f = (/i,/2) takes on values in a triangular 
domain with the three collinear states in its corners. 

What kind of bond order might one expect in the in- 
termediate phase? To answer this question, let us again 
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FIG. 8: The distribution of bond vector of the two sublattices 
fA and fs in the nematic phase. The simulated system has 
(a) 8 layers and (b) 6 layers of tetrahedra in one sublattice. 
The bond vector f has been normalized such that the three 
collinear states, blue, red, and green, are at vertices ( — 1,0), 
(i, ^), and (i, respectively. 



use the equivalence between a ferromagnetic Ji and an 
antiferromagnetic J3. The latter promotes antiparallel 
orientations for spins 3 and 3' (Fig. [IJ, which means — 
for a collinear state of spins — that one of the bonds 03 
and 03' is frustrated and the other is satisfied. (Bergman 
et al£^ showed that such states — satisfying the "bend- 
ing rule" for frustrated bonds in zero applied field — are 
also favored by quantum fluctuations of spins.) In other 
words, adjacent tetrahedra will be in states of different 
color. This is reminiscent of the antiferromagnetic Potts 
model with 3 states: red, green and blue in Fig. [7] A 
collinear state of the pyrochlore antiferromagnet is fully 
specified by the global spin director and the colors of 
all tetrahedra. Note however that colors of tetrahedra 
are not completely independent: the number of satisfied 
bonds (S;-Sj = — S 2 ) must be even along any closed loop. 
Nonetheless, the parameterization in terms of Potts vari- 
ables serves a useful purpose. One of the phases of the 
antiferromagnetic Potts model on a bipartite lattice has 
a broken sublattice symmetry (BSS): one sublattice is 
dominated by one color, while the other is randomly pop- 
ulated by the two remaining colors i 25 ' 26 With this state 
in mind, we have measured the average bond variables in 
the intermediate phase in the Monte Carlo simulations. 
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FIG. 9: Histogram of 17 distinct collinear layered structures 
obtained by replica-exchange Monte Carlo simulation. The 
system has 16 x 4 3 spins. The configuration number labels 
17 topologically distinct layered states subject to the periodic 
boundary condition. 



The Monte Carlo averages of the bond doublet ^ for 
sublattices A and B are shown in Fig. [5) The value of 
f for sublattice A is narrowly distributed in the vicinity 
of the collinear blue state, indicating that all tetrahe- 
dra of sublattice A are in this state. There are no blue 
tetrahedra on sublattice B, as one might expect from the 
analogy with the antiferromagnetic Potts model. For the 
BSS phase, where each site is red or green with equal 
probabilities, one expects a continuous distribution of f 
in the middle of the opposing edge of the triangle con- 
necting the green and red corners. Instead, we find that 
sublattice B has discrete fractions of red tetrahedra: e.g. 
0, 1/4, 1/2, 3/4, and 1 in a system with 8 layers of tetra- 
hedra in one sublattice [Fig. [UJa)]. 

This discreteness is a finite-size effect: an examina- 
tion of individual microstates shows that the intermedi- 
ate phase has a layered structure for bond variables on 
sublattice B: tetrahedra within the same layer in the 
xy plane have the same color. The origin of the layered 
structure on one of the sublattices can be traced to the 
same constraint on the colors around a closed loop. See 
Appendix |B] for details. For example, the simulated sys- 
tem of Fig. [H^a) contained 8 layers of tetrahedra within 
a sublattice. If the layers could be colored red and green 
independently of one another, one would expect to find 
the fractions of either color proportional to 1/8. How- 
ever, periodic boundary conditions create constraints on 
the number of satisfied bonds in the direction perpendic- 
ular to the layers, so that each lattice can only have an 
even number of layers of either color. Hence the fractions 
proportional to 1/4. Similarly, for a system in which each 
sublattice has 6 layers of tetrahedra, the fraction of red 
layers is 0, 1/3, 2/3, and 1 [Fig. Mb)). 

To verify this observation more directly, we performed 
a replica-exchange Monte Carlo simulation on a system 
with 4 3 conventional cubic cells. 16 x 4 3 spins are divided 
into 8 layers of tetrahedra in each sublattice. A partic- 
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ular layered state with collinear spins is described by a 
sequence of Ising variables {ci, cx 2 , • • • ^s} ( see Appendix 
ICj) . With periodic boundary conditions, 17 distinct con- 
figurations are used in a replica-exchange Monte Carlo 
simulation. The Ising sequences corresponding to these 
17 layered states are listed in Table I. In each exchange 
cycle, a fixed number of Metropolis sweeps are performed 
on individual replicas of the system, each of which corre- 
sponds to a particular layered state. Then different repli- 
cas are exchanged according to detailed balance, thus en- 
suring thermodynamic equilibrium. A histogram of the 
occurrence of the 17 configurations in a chosen replica is 
shown in Fig. [9j The almost equal probability of occur- 
rence implies a vanishing spin order after averaging over 
the different configurations. 

The layered structure of the intermediate phase spon- 
taneously breaks the rotational and translational sym- 
metries of the pyrochlore lattice. A collinear Neel order 
exists within an individual layer of tetrahedra but not 
across the layers if the colors on one sublattice are in- 
deed random. At the mean-field level, the collinear states 
in the partially ordered phase belong to a larger class of 
(generally non-collinear) layered states with the same ex- 
change energy. A discussion of the general layered states 
is presented in Appendix [Cj As already mentioned pre- 
viously, since collinear spins tend to have softer magnon 
spectrum, those layered states with collinear spins are 
favored by thermal fluctuations. 

The two phase boundaries enclosing the intermediate 
phase are both discontinuous transitions. The critical 
temperatures determined by the mixed-phase method 17 
arc linear in J 2 : T cl ~ 1.87| J 2 |S 2 and T c2 ~ 0.26| J 2 \S 2 
as T — > 0. Our numerical simulations seem to indicate 
that the intermediate phase is globally stable in the tem- 
perature regime T c2 < T < T c \: in the mixed state, 
the collinear phase gradually takes over the entire lat- 
tice. We do not have analytical arguments to back up 
the global stability of the intermediate collinear phase: 
such an analysis would require knowledge of the free en- 
ergy of the magnetically ordered phase, which has not 
yet been obtained. 



IV. LOCAL STABILITY OF THE PARTIALLY 
ORDERED PHASE 

Even an analysis of the local stability of the partially 
ordered collinear phase is not exactly straightforward. 
The standard large-5* method of computing the magnon 
contribution to the free energy fails because of the exis- 
tence of unstable modes with a negative stiffness at zero 
temperature. The instability merely reflects the fact that 
the collinear states are not a local minimum of energy |T]). 
The instability is avoided at a (sufficiently high) finite 
temperature: the free energy of spin fluctuations con- 
tributes a positive term to the spin stiffness. In this Sec- 
tion we analyze the local stability of the collinear phase. 



(a) (b) 



h 



@ 




® 




») 


<§> 




C 


*) 


<@> 




<J 


<§> 
















») 


( ? ) 




( 


») 


<§> 















(h, h, 0) / 

FIG. 10: (a) The energy dispersion of the spin-wave band 
with unstable modes, (b) Regions in momentum space q = 
2ir(h, h, I) where the spectrum of energy fluctuations has neg- 
ative eigenvalues Aq. Ji = 1, J2 ~ —0.1. 



A. Unstable modes 

To analyze the stability of a collinear state, we ex- 
press the energy of the system in terms of transverse 
spin fluctuations er^ _L n. By substituting S, rs 5*(1 — 
erf /2S 2 )r]ih + <Ji into Eq. (JXJ) we obtain a spin-wave 
Hamiltonian in the harmonic approximation, 

H i2) =E L + (Jx - 2J 2 ) °f + \ J v ai ■ ( 6 ) 

i i,j 

where El is the energy of the layered state. The Ising 
variables {7^} specifying the direction of a spin are absent 
from the harmonic Hamiltonian ([6]). They affect the dy- 
namics of the system through the canonical commutation 
relations for the transverse components of the spins. 

The quadratic form © must be positive definite to 
guarantee stability of the collinear state. Its eigenvalues 
A are obtained by making the Fourier transform and then 
diagonalizing a 4 x 4 matrix (the pyrochlore lattice is an 
fee with a basis of 4 sites) : 

A^ = (J x - 2J 2 ) + A°, (7) 

where Aq are eigenvalues of J mn (q) defined in Sec. Ill Bl 
The dispersion has degenerate zero modes along lines 
q = 27r{l,/i, 0} corresponding to magnetic spirals along 
one of the three cubic axes. These spirals belong to the 
degenerate manifold of non-collinear layered states dis- 
cussed in Appendix [C] Furthermore, there are regions 
in momentum space with A q < 0, as shown in Fig. 
[T0l The most unstable modes are found at wavevectors 
q* = 2n{h*,h*,0} with h* given by Eq. g]). For small 
J2/J1, the lowest eigenvalue is 

^ = (28 - 16V3) ^ + ^(56V3 - 97) f J -ff + ■ ■ ■ . 

Since A m j n < for a ferromagnetic J 2 , the collinear 
ground states are unstable at zero temperature. 
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FIG. 11: (a) energy density e and (b) nematic order parameter 
as a function of temperature obtained using Monte Carlo sim- 
ulations and a Hartree-Fock self-consistent calculation. The 
calculation was done with J2 = — O.OUi. The dashed line 
is a linear fit to the Monte Carlo data. Note that the tran- 
sition temperature obtained from Monte Carlo simulation is 
T c2 « 0.076 \J 2 \S 2 . 



B. Hartree-Fock calculation 

At finite temperatures the collinear layered states are 
stabilized by thermal fluctuations. To demonstrate this, 
we go beyond the harmonic term of the classical Holstein- 
Primakoff expansion and consider the interactions be- 



tween spin waves 



27 



1 

8S 2 



1 



ai ■ a, (<r 2 + a 2 .) .(8) 



Since the system is unstable at the harmonic order, a per- 
turbation expansion based on the quadratic Hamiltonian 
© is not possible. Instead, following Hizi and Henley, 2 ' 
we construct an effective (mean-field) quadratic Hamil- 
tonian 



n 



MF 



(9) 



that provides the best approximation to TL^ +H^. To 
this end, we use the standard mean-field recipe to decou- 
ple the quartic Hamiltonian. We first write every possible 
pair of operators in 7v 4 ) in terms of its thermal average 
plus a fluctuation term. Dropping terms quartic in the 
fluctuations yields the quadratic form © with the fol- 

~ (2") 

lowing coefficients : 
( Jl - 2J 2 ) + J2k J ik (Vi Vk G k k - G ik ) (i = j), 



2 J y 



1 + gW* Vj G i 



2s 1 



(G H + G js ) (i^j).(10) 
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FIG. 12: Stability boundary T* obtained using the Hartree- 
Fock calculation and the Monte Carlo simulations. The error 
bars shown for the Monte Carlo data are equal to the tem- 
perature step AT used in the simulation. 



Here Gij = 



is the correlation function of 



spin fluctuations calculated self-consistently in the ther- 
mal ensemble of the mean-field Hamiltonian (fH 



Gij — 



I Da afa* e 



J Da, 



-/?Hmf 



(ii) 



Numerically, an iteration process is used to obtain 
the correlation functions Gij. After self-consistency is 
reached, the energy of the magnet is given by 

Emf = E L + 2 / ](Ji — 2 J 2 ) Gu + )] Jjj G^ (12) 

i i,j 

+ t^t ^ IpiVji^uGjj + Gij) - G i j{G li + Gjj) 

Fig. [11] (a) shows the computed energy density as a func- 
tion of temperature. The result agrees very well with that 
obtained from Monte Carlo simulations. Both the simu- 
lation and calculation were done for J2 = —0.01 J% on a 
pyrochlore lattice with a size of 16 x 4 3 spins and peri- 
odic boundary condition on each side. The self-consistent 
method can also be used to compute the nematic order 
parameter. For n = +z, the tensor (S^S^) becomes 
diagonal with elements (S X S X ) — (S y S y ) — 2G and 
(S Z S Z ) =1-2(5, where G = J2i G n/ N - The nematic 
order parameter is then 



Q 



NS 2 



5> 



(13) 



The result is shown in Fig. [TT] (b) and the agreement 
with that obtained from Monte Carlo simulation seems 
satisfactory: the discrepancy between the two methods 
is less than 3%. The nearly saturated nematic order 
parameter Q observed in Monte Carlo simulations im- 
plies er 2 <C 1, justifying the Holstein-Primakoff expansion 
about the collinear state. 
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FIG. 13: Renormalized nearest neighbor bonds of the red- 
and-green state in the mean-field calculation. The renormal- 
ized first-neighbor exchange constants: Ji — K\ — K2 (dashed 
bonds), Ji — K\ (dash-dotted bonds), and Ji — K2 (solid 
bonds). 



Below a certain temperature T* the energy spectrum 
of spin waves acquires some negative eigenvalues and the 
collinear phase gives way to the low-temperature ordered 
state. Since the transition is first order, the E — T dia- 
gram exhibits hysteresis. The thermodynamic transition 
takes place at a temperature T c2 > T*, at which the 
collinear phase is still locally stable. 

The dependence of T* /\J%\ on the ratio \J 2 \/ Ji ob- 
tained from the Hartree-Fock calculation is shown in Fig. 
1121 The points collapse perfectly on a linear curve imply- 
ing a scaling relation T* ~ J 2 / J\ . A numerical estimate 
of the stability boundary T*(J 2 ), obtained as the lowest 
temperature at which the intermediate phase was still 
observed in Monte Carlo runs, is also plotted in Fig. [T2l 
the result is in satisfactory agreement with that of the 
mean-field calculation. 



8J01 


= SJ23 = 


-K x - K, 


SJ02 


= SJ 31 = 




SJ03 


= 5J 12 = 


-K 2 



number of variational parameters reduces to 2. This 
is so because zero modes satisfy Ei=o Gi = 0, hence 
(<7o<7i) = — ((Jq) — (0002) — (0003). It follows then that 
SJqi = 5J02 + SJ03. We parameterize the exchange renor- 
malizations in terms of K\ and K 2 such that 



(15) 



on the red sublattice. 

We then compute the spectrum and the eigenmodes 
of energy fluctuations with the renormalized exchange 
interaction. The two zero-energy bands that were flat in 
the absence of J 2 and K t now acquire a dispersion; one 
becomes gapped (A^ is strictly positive), while the other 
has a vanishing energy at the wavevector qo = 27r(0, 0, 1). 
This zero mode corresponds to a global rotation of spins. 
Correlation functions are dominated by fluctuations in 
the lowest band in the vicinity of qo- For small k, the 
energy eigenvalue is 



A^+k^i^fci 



K 2 )k 2 z 



(16) 



where k\_ 



In order to obtain the correlations Gij , we need first to 
obtain the eigenmodes. To this end, we use an orthonor- 
mal basis of the two zero modes of J\ for given values of 
k. We then treat Ki and J 2 as perturbations and use de- 
generate perturbation theory to obtain the eigenmodes. 
To the lowest order in k, they are 



C. Analytic results: red-and-green state 

An analytical derivation of the stability temperature 
T* ~ J 2 /Ji is difficult to obtain for the most general 
layered state. We have evaluated the stability for the 
simplest state of this kind, where all of the layers have 
the same colors. A state of this sort (sublattice A is 
red and sublattice B is green) was studied in Ref. [H. 
This particular state has a higher symmetry than a typ- 
ical layered structure: the color variables violate only 
the inversion symmetry exchanging the two sublattices 
of tetrahedra. 

In the mean-field Hamiltonian (jTUJ) , the main effects of 
the quartic interaction TL^ is to renormalize the first- 
neighbor exchange J\ 
bond-dependent : 



to Jij = Ji 



SJij, which is now 



G 



30 



2T)iT)jGij). 



(14) 



Assuming that exchange renormalizations SJij respect 
the symmetries of the red-and-green state, we have 3 in- 
dependent variational parameters SJ01, SJo 2 , and <5Jo3 
(Fig. [TB"]) . If we further assume that the correlations 
Gij are dominated by the pyrochlore zero modes, the 



w (q +k) = -i/2 - (k x - k y + k z )/16, 

wi(q + k) = +1/2 - i(k x + k y + k z )/16, 

w 2 (qo + k) = -l/2-i(k x + k y -k z )/16, 

u 3 (q Q +k) = +i/2-(k x -k y -k z )/m. (17) 

As can be easily checked, the total spin of a tetrahedron 
E™ o>n = E™ Wme l(qo+k) r " 1 = at this order of k. The 
spin correlation function is 



G -» = ^7 E A7 <re(q) W " (q) el(qHr " 
q 



-r„) 



(18) 



where N 1 — N/4 is the number of unit cells, and m, n are 
sublattice indices. By expanding to the second order of 
k and using (fTl)) . we obtain the following self-consistency 
equations for K\ and K 2 



AN'S 2 ^ 2K X k 2 , 

k - 1 - 

JiT 

2N'S 2 ^ 2Ki k'i + (8\J 2 \ + K 2 )k 2 



E 



(8\J 2 \+K 2 )k 2 
k 2 ± 



(19) 
(20) 



Although these equations can be solved numerically, 
we are interested in an approximate solution of K\ and 
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FIG. 14: Variation of spin- wave energy A in unit of | J%\ along 
the q = 2n(h, h, 1) line. The calculation was done with a Ji = 
— O.OUi. The curves correspond to temperatures T l\Ji\ = 
0.018, 0.0165, 0.01526, 0.0145, 0.013 (from top to bottom) 
T* = 0.01526| J2 15 2 corresponding to the temperature where 
the q = mode becomes unstable. 



the instability, whereas the proper calculation at a finite 
temperature should involve the free energy. We do this 
next. 



D. Stability boundary: red-and-green state 

We now provide an estimate of the stability temper- 
ature T* by computing the magnon contribution to the 
system free energy An expression (|D9[) for the change of 
free energy associated with an unstable mode is derived 
in Appendix [D] Here we apply the result to the red- 
and-green state. We consider the most dangerous modes, 
namely those with wavevectors near q* = 2n{h*, h* , 0} 
where h* rs 3/4. In the presence of such an unstable 
mode with amplitude <f> superimposed on the red-and- 
green state, the free energy changes by an amount given 
by 



AF 



(a* 



S 2 



(23) 



Ki in the low-temperature regime, T -C | J2IS 2 . Since the 
effective spin stiffness K is generated by thermal fluctua- 
tions, they are expected to be small compared to J2. To 
the lowest order we neglect K\ and K2 in Eq. (fTi?)) and 
obtain 



JiT 



32|J 2 |S 2 ' 



(21) 



On the other hand, because the integral for K2 is diver- 
gent as ifi -> 0, we must keep K\ in Eq. (f20"|) . Substi- 
tuting the result for K\ into Eq. (|20|). we obtain 



K 2 



3y/2S 



(22) 



to the lowest order in T. 

These results provide a glimpse into the physics of the 
transition between the intermediate and low-temperature 
phases. Fig. RBI shows the rcnormalized dispersion A qo+ k 
(fT6|) along the line q = 2w(h, h, 1) at various tempera- 
tures. As the temperature decreases, a dip of the disper- 
sion curve starts to develop at h » 0.2. Eventually this 
local minimum touches zero at the critical temperature 
T C 2\ below T* the collinear state is unstable: it decays 
by emitting spin waves with q w 27r(l/4, 1/4, 1), which 
is related to 27r(3/4, 3/4, 0) by a reciprocal lattice vector. 

It should be noted that the scenario displayed in Fig. 
1141 is only a qualitative description of the real transi- 
tion. Our self-consistent treatment only takes into ac- 
count spin waves close to the qo = 27r(0,0, 1) Goldstonc 
mode. This is valid at temperatures well above T* since 
these spin waves are the lowest-energy excitations of the 
magnet. However, as T — > T C 2, spin waves with wavevec- 
tors q « 27r(3/4, 3/4, 0) become soft and should also be 
included in a self-consistent calculation. Additionally, 
we have studied the energy of spin waves as a proxy for 



where the correlation function G mn is given by fTB]). 
and A mn is the perturbation to the mean-held Hamil- 
tonian 7Ymf caused by the unstable mode. In our case, 
the real-space eigenvector of the unstable mode with 
q* = 2n(h*,h*,0) is 



m„(r) = U* [xcos(q* ■ r) + y sin(q* • r)] , 



(24) 



where the corresponding momentum-space eigenvector 
for q* is 

U* = (cos9, - sin 9, - sin 9, cos6>)/\/2, (25) 

with 8 w 0.27-7T and weakly dependent on J 2 . We write 
the energy of the unstable mode as A* = —7 \J%\, where 
7 w 0.2 is a dimensionless number. The change of free 
energy is then 



AF/0 2 = - 7 | J 2 \S : 



47V' ^ 2K1 k\ + (8\J 2 \ + K 2 )k 2 z 



where 



^ A m „u* (k)w m (k)e 



i(qo+k)-(r„ 



(26) 



(27) 



Since in most cases Ak ~ A + 0(k 2 ) for h* = 1/4, we 
neglect the k dependence of Ak in the following as a 
lowest order approximation. With the aid of Eqs. QT 
and (|20[) . the integral evaluates to 



V2A 
16tt 



(28) 



The condition AF = thus gives an estimate of the 
stability temperature 



V y^A ) Ti 



(29) 
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This expression overestimates (by a factor of about 10) 
the stability temperature compared with numerical re- 
sults. However, as mentioned previously, the discrepancy 
is due to the fact that we neglect contributions from the 
unstable modes themselves when approaching the tran- 
sition temperature. Those modes with wavevector cen- 
tered about the 12 unstable q* = 2n{h*, h* , 0} become 
extremely soft as T — > T* and should be included in 
the calculation in a self-consistent way. Nevertheless, 
provides an upper bound of the stability boundary 
and gives a scaling relation consistent with the numerical 
data. 



V. DISCUSSION 

We have studied the classical Heisenberg antiferro- 
magnet on the pyrochlore lattice with first and second- 
neighbor exchange interactions. Ferromagnetic second- 
neighbor exchange J 2 < is frustrated and lifts the vast 
degeneracy of the nearest-neighbor model only partially, 
setting stage for a nontrivial phase diagram in the ( J2, T) 
plane. We have used a combination of Monte Carlo sim- 
ulations and analytical calculations to characterize the 
phases of this model. In our opinion, the low-temperature 
phase, discussed previously by Tsuneishi et al.^ is the 
incommensurate, and likely non-collinear, ordered phase 
predicted earlier by Reimers et alJ£ A full characteriza- 
tion of its magnetic order remains to be done, and its 
fate in the presence of strong quantum fluctuations is an 
interesting topic for future study. 

Our simulations have uncovered the existence of an- 
other, partially-ordered phase at intermediate tempera- 
tures for a weak enough \ J%\. In the intermediate phase, 
the spins are on average collinear, which is manifested by 
a nonzero nematic order parameter. The order is fully 
characterized by a combination of a global nematic di- 
rector n and a 3-state Potts variable (color) on every 
tetrahedron indicating the location of frustrated bonds 
(Fig. [7]). The second- neighbor interaction J2 < acts 
like an antiferromagnetic Potts coupling forcing unlike 
colors on neighboring tetrahedra. 

The color structure of this phase resembles the ordered 
state with broken sublattice symmetry (BSS) of the an- 
tiferromagnetic Potts models one sublattice of tetra- 
hedra is dominated by one color (say, blue) while the 
other exhibits a mixture of the remaining two colors (red 
and green). However, unlike in the BSS state, the two 
colors on the second sublattice are not distributed in a 
completely random way: they form uniform layers in the 
plane associated with the colors (in this case, xy). The 
colors of individual layers appear to be random, hence 
partial order. 

The partial order can be described by an individual Z2 
variable <7j for each such layer — in addition to a global di- 
rection of the spins n and the color of the other sublattice. 
States with different sets of {ui\ are local minima of the 
free energy. Accessing one such minimum from another 



by means of a uniform rotation of spins within one layer 
of tetrahedra requires climbing over a free-energy barrier 
that grows as the number of spins in that layer and thus 
becomes impossible in the thermodynamic limit. A more 
plausible route to changing the color of a layer is by nu- 
cleating a bubble of the opposite <7j, which will grow if 
the new state has a lower free energy once the bubble is 
large enough for the gain in bulk energy to outweigh the 
cost in interface energy. Since the distinct layered states 
are not related by symmetry, their free energies are gen- 
erally different and the nucleation route may well lead to 
a selection within this set of states. Since such nucleation 
can go along with large energy barriers, it can be tricky 
to observ o 28 ' 29 , and indeed we have not found it in our 
simulations. 

It is worth stressing that the ideal collinear states 
do not minimize the exchange energy — either globally 
or locally. They owe their stability to thermal fluctu- 
ations, which effectively renormalize exchange couplings 
and turn these spin configurations into minima of the free 
energy. As the temperature falls, the couplings return to 
their bare values and the collinear states become locally 
unstable at a temperature T* = C(j|/Ji), in agreement 
with our Monte Carlo simulations. The most unstable 
spin-wave mode has approximately the same wavenum- 
ber as the low-temperature incommensurate magnetic or- 
der. The simulated phase transition is strongly discon- 
tinuous. 

Simulations on the high-temperature side show that 
the intermediate phase persists up to a temperature 
0(J2)- A discontinuous phase transition takes it into 
the paramagnetic phase. The presence of strong local 
spin correlations in the paramagnetic phase means that 
the effect of third-neighbor couplings J3 (but not of J3, 
see Fig. [T]) is equivalent, up to a change of sign, to that 
of the second-neighbor coupling, at least to the first or- 
der. Therefore we expect that the state of our system 
depends on these couplings mostly through their differ- 
ence J 2 — Jz- If correct, this observation would extend 
the results of our study to a broader class of pyrochlore 
antiferromagnet with both J2 and J3 present. 
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FIG. 15: A fragment of a collinear state. Frustrated bonds 
are shown as colored dashed lines. 



APPENDIX A: EQUIVALENCE OF J 2 AND -J 3 
IN THE STRONGLY CORRELATED 
PARAMAGNET 

At temperatures T <C J\S 2 spins on every tetrahedron 
approximately satisfy the constraint 

£> = 0. (Al) 

i=0 

Consider the effective magnetic field on the site labeled 
3' in Fig. [[J 

H 3 , = -dH/dS 3 , = - J! S - J 2 (Si + S 2 ) + . . . (A2) 

where we have explicitly written out the contributions 
from the spins of the adjacent tetrahedron 0123. Let us 
now turn off the second-neighbor exchange, J 2 = 0, and 
turn on the third-neighbor coupling J 3 (Fig. [T]). Doing 
so changes the effective field to 

- JiS - J 3 S 3 « -(Ji - J 3 )S + J 3 (Si + S 2 ), (A3) 

where we used the constraint (jAip . In this setting, a 
comparison of Eqs. (|A2[) and (|A3|) shows that adding a 
third-neighbor coupling J3 is indeed energetically equiv- 
alent to the second-neighbor exchange of the same mag- 
nitude and opposite sign. This result does not extend to 
the excited states, which violate (jAll) . so that the physics 
of fluctuations need not be simply related. 

APPENDIX B: CONSTRAINT ON COLORS 
(BOND VARIABLES) 

Consider the hexagonal loop abcdef shown in Fig. [T51 
Suppose that tetrahedra of sublattice A are in the blue 



state and that one tetrahedron B\ of the other sublattice 
is red. Then it can be seen that tetrahedron i? 2 , which 
has the same z coordinate, must also be red. This can 
be proved as follows. In the collinear state, spins can be 
represented by a Ising variable, i.e. Si = Srjiii and n is 
an arbitrary unit vector. Obviously, the product of the 
six bond variables rjiijj on the hexagon loop is +1, i.e. 

{VaVb)(VbVc) ■ ■ ■ {VeVf)(VfVa) = +1. (Bl) 

Among the six bonds, sublattice A contributes two an- 
tiferromagnetic and one ferromagnetic bond; this makes 
its total contribution +1. Therefore the product of the 
three bonds on sublattice B must be +1 as well. We 
know that rjfrj a = +1 {B\ is red) and rjbrjc = — 1 (_B 3 is 
not blue) . Hence r/^r/ e = — 1 , which means that _B 2 is not 
green. Since -B 2 is not blue, it must be red. 

Thus, if sublattice A is blue, the above proof shows 
that sublattice B has uniform color green or red in each 
layer z — const. However, the color of individual B layers 
are random. This is similar to the BSS phase of antifer- 
romagnetic Potts model where individual sites on one 
sublattice have random colors. 



APPENDIX C: MAGNETIC STRUCTURE OF 
THE LAYERED STATE 

At the mean-field level, the collinear BSS-like states of 
the intermediate phase are degenerate with a larger class 
of layered state with non-collinear spins in general. Here 
we describe the magnetic structure of the general layered 
state. 

At temperatures well below the Curie- Weiss constant, 
the magnetic state of a tetrahedron is determined by 
three staggered magnetizations L,-, where Li = (So + 
Si — S 2 — S 3 )/4, and so onJ^ Here we choose to specify 
the Neel vectors of layers belonging to sublattice A. Be- 
cause each spin is shared by two tetrahedra from different 
sublattices, the magnetic state of tetrahedra of sublattice 
B is encoded in the staggered magnetizations of the four 
surrounding tetrahedra of sublattice A. 

Choosing the normal of the layers to be the z axis, the 
staggered magnetizations of a tetrahedron with inplane 
coordinate vector in the A:-th layer are 

L 1= L 2 =0, L 3 = Sn k e^ r \ (CI) 

where = 27r(l,0) or 27r(0,l) which are equivalent 
with respect to the two-dimensional (2D) square lattice 
of tetrahedra with the same z coordinate, and is an 
arbitrary unit vector. Fig. 1161 shows an example of the 
2D Neel order. Once the spin order within the layers 
is specified, the magnetic structure of a general layered 
state is described by a sequence of the unit vectors . 

The bond order of the layered state is similar to the 
BSS state of three-state Potts model. In the example 
given above, sublattice A is in the collinear blue state 
while tetrahedra in sublattice B in general have coplanar 



13 




FIG. 16: Non-collinear layered state projected along the x- 
axis. The normal of the layers are parallel to the z axis. 
Frustrated bonds are shown as dashed lines. 



spins; their bond order is determined by the Neel vectors 
of the two A layers enclosing it: 

45 2 /l V3\ „ \ 

fs = 7fl2'— nfc ' nfc+1 J- (C2) 

For arbitrary n^, the bond vector is spreads uniformly 
on the edge of the triangle domain which connects the 
two vertices corresponding to the red and green states; 
the average color is again yellow. 

The energy of a layered state is independent of the 
direction rij. of spins in the individual layers: 

E L = -N(Jx - 2 J 2 )S 2 . (C3) 

This energy corresponds to the extrapolated zero- 
temperature energy density sl = — | 1 / 3 of the nematic 
phase in Fig. O Although all layered states are degen- 
erate at the mean- field level, thermal fluctuations appar- 
ently prefer the collinear ones as shown by the Monte 
Carlo simulations. 

In the collinear layered states, a common direction n 
is selected and the unit vectors hk — > Uk n with the Ising 
variable = ±1. The magnetic structure of a layered 
state is then specified by a sequence of Ising variables: 
{ci, 02, • • • , (Tfe, • • •}. For a pyrochlore lattice with 8 lay- 
ers of tetrahedra in each direction, there arc 17 distinct 
layered states that are not related to each other by trans- 
lations and inversions of the Ising variables (Table I). 
However, some of these states may be related by other 
symmetries of the lattice. For example, both states 1 and 
17 represent the red-and-green state. 



No. 


Ising sequence 


9 


M h + + 


1 


++++++++ 


10 


-+-+-+++ 


2 


-+++++++ 


11 


- + H h H h 


3 


h + + + + + 


12 


- h + + + 


4 


— 1 h + + + + 


13 


1 h + + 


5 


-++-++++ 


14 


- - - + + -+ + 


6 


- + + H h+ + 


15 


1 h - + + 


7 


+ + + + + 


16 


1 h H h 


8 


1 h + + + 


17 


- H 1 1 h 



TABLE I: Ising sequences of the 17 distinct layered states for 
a pyrochlore lattice with 8 layers subject to periodic boundary 
conditions. 



APPENDIX D: FREE ENERGY OF THE 
UNSTABLE MODE 

Below we derive a Holstein-Primakoff Hamiltonian for 
spin waves in a state with nearly collinear spins and com- 
pute the free energy to test the local stability of the 
collinear state. We focus specifically on the previously 
identified unstable modes. 

In a general state with noncollinear spins Si = Shi, 
we introduce a local reference frame defined by three or- 
thonormal vectors: ef, ef , and n^. A small deviation 
from this state can then be expressed using the Holstein- 
Primakoff expansion: 

S< = S(l - |J) fii + E a? ef + 0(a 3 ), (Dl) 

ct—x.y 

Here <Ji = (of, of) whose components denote fluctua- 
tion along the two orthogonal local axes. The exchange 
Hamiltonian then becomes 

h = J2 Ja a* ■ % + E E H f < ■ (° 2 ) 

ij id a, (3 

with 

Hf = \7\^ a J ^- hkda " f.Z^ (D3) 

For collinear states in the nematic phase, = ^z, and 
ef = x, e| = y, Eq. (|D3|) reproduces the harmonic 
Hamiltonian ([6]). 

Next we compute the increase in the free energy result- 
ing from a small deviation from a collinear state in the 
direction of an unstable mode of the bare Hamiltonian 
<j6j> . Let 4> be the amplitude of the unstable mode and 
{mi} the corresponding (normalized) real-space eigen- 
vector. The deformed spin configuration is 

Si = Sh % = Sirn(l - 2 m2/2) + S^m, (D4) 

Given this local spin axis , there is arbitrariness in the 
choice of the other two unit vectors ef . In order to apply 
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the perturbation method, we choose: 



(D7) 



:«U 



2 

: (mf) 2 



- ZT)i <j> m^.(D5) 



Substituting these expressions into Eq. (|D3[) . we obtain 
H°f = 5 a(} H[f + <j) 2 A?f, where H^f is the magnon 
Hamiltonian of the collinear state and the perturbation 
i s given by 

A «/3 = f ~\ Efe [ m * ' m fc ~ kViVki-m-i + m 2 k )]5 a P, (i 

tJ \ \ JijiViVj m f m j — ^( m ? m i + m j' m j)]> (* 

; ~ (D6) 

Since the bare harmonic Hamiltonian contains unstable 
modes as discussed in Section ITVl we replace by the 
one renormalized by spin-wave interactions, given in Eq. 
(HHJ). We may then approximate the free energy of the 
system as 



-P F » p -f3{E L +A'S 2 <j> 2 ) 

n -/3£ [ff <2 W°-,+^ 2 A Q V a <7' 5 



Here is energy of the layered state, Z is the parti- 
tion function of the renormalized Hamiltonian Ti^ , (• • •) 
means Boltzmann averaging with respect to the Hamilto- 
nian 7t(( 2 ), and A* < is the bare energy of the unstable 
mode nii. The prime in the integral indicates that we 
only integrate out the low-energy magnons close to the 
Goldstone mode of the collinear state. 

Upon expanding the fluctuations in terms of spin- wave 
eigenvectors erf = J2 n £n u n,i, we obtain the spin corre- 
lation: 
= 3) 



+ 3) 



cr,- a 



3 ' 



— () a P n 

n,j — u ^13 ■ 



(D8) 



Here A n is the energy of the n-th eigenmodc of the renor- 
malized Hamiltonian . Substituting this result back 
into Eq. (|D7|) yields the free energy (|23|) associated with 
the unstable mode <j>, 

F « const + (k*S 2 + Y,G^3 A Ji)^ + 0(^ 4 ), (D9) 



where Ay = Y, a 



a a 
ij ' 



1 R. Moessner and A. P. Ramirez, Physics Today 59 (2), 24 
(2006). 

2 G. H. Wannier, Phys. Rev. 79, 357 (1950). 

3 R. M. F. Houtappel, Physica 16, 425 (1950). 

4 S. T. Bramwell and M. J. P. Gingras, Science 294, 1495 
(2001). 

5 R. Moessner and J. T. Chalker, Phys. Rev. B 58, 12049 
(1998). 

6 S. V. Isakov, K. Gregor, R. Moessner, and S. L. Sondhi, 
Phys. Rev. Lett. 93, 167204 (2004). 

7 C. L. Henley, Phys. Rev. B 71, 014424 (2005). 

8 C. Kittel, Phys. Rev. 120, 335 (1960). 

9 O. Tchernyshyov, R. Moessner, and S. L. Sondhi, Phys. 
Rev. B 66^ 064403 (2002). 

S.-H. Lee, C. Broholm, T. H. Kim, W. Rat cliff II, and 
S.-W. Cheong, Phys. Rev. Lett. 84, 3718 (2000). 

1 J.-H. Chung, M. Matsuda, S.-H. Lee, K. Kakurai, H. Ueda, 
T. J. Sato, H. Takagi, K.-P. Hong, and S. Park, Phys. Rev. 
Lett. 95, 247204 (2005). 

2 J. N. Reimers, A. J. Berlinsky, and A.-C. Shi, Phys. Rev. 
B 43, 865 (1991). 

3 D. Tsuneishi, M. Ioki, and H. Kawamura, J. Phys. Con- 
dens. Matter, 19, 145273 (2007). 

4 O. Tchernyshyov, O. A. Starykh, R. Moessner, and A. G. 
Abanov, Phys. Rev. B 68, 144422 (2003). 

5 C. L. Henley, Phys. Rev. Lett. 62, 2056 (1989). 

6 C. Pinettes, B. Canals, and C. Lacroix, Phys. Rev. B 66, 



024422 (2002). 

17 M. Creutz, L. Jacobs, and C. Rebbi, Phys. Rev. D 20, 
1915 (1979). 

18 G.-W. Chern, C. Fennie, and O. Tchernyshyov, Phys. Rev. 
B 74, 060405(R) (2006). 

19 A. N. Yaresko, Phys. Rev. B 77, 115106 (2008). 

20 C. L. Henley, Phys. Rev. Lett. 62, 2056 (1989). 

21 P. M. Chaikin and T. C. Lubensky, Principles of condensed 
matter physics (Cambridge University Press, Cambridge, 
2000). 

22 K. Hukushima and K. Nemoto, J. Phys. Soc. Jpn. 65, 1604 
(1996). 

23 S. Trebst, D. Huse, and M. Troyer, Phys. Rev. E 70, 
046701 (2004). 

24 D. L. Bergman, R. Shindou, G. A. Fiete, and L. Balents, 
Phys. Rev. B 74, 134409 (2006). 

25 G. S. Grest and J. R. Banavar, Phys. Rev. Lett. 46, 1458 
(1981). 

26 S. Lapinskas and A. Rosengren, Phys. Rev. Lett. 81, 1302 
(1998). 

27 U. Hizi and C. L. Henley, J. Phys. Cond. Matt. 19, 145268 
(2007). 

28 A. J. Leggett, Phys. Rev. Lett. 53, 1096 (1984). 

29 P. Schiffer, M. T. O'Keefe, M. D. Hildreth, Hiroshi 
Fukuyama, and D. D. Osheroff, Phys. Rev. Lett. 69, 120 
(1992). 



